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Abstract 



We have calculated inclusive two-jet cross sections in next-to-leading order QCD for 
^ , direct photoproduction in low ep collisions at HERA. Infrared and collinear singular- 

?H ' ities in real and virtual contributions are cancelled with the phase space slicing method. 



Analytical formulas for the different contributions giving the dependence on the slicing 
parameter are presented. Various one- and two-jet distributions have been computed 
demonstrating the flexibility of the method. 
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1 Introduction 

Electron-proton scattering at HERA proceeds dominantly through the exchange of photons 
with very small virtuality {Q^ ~ 0). An important fraction of the 7p collisions contains high 
transverse momentum {Ex) processes. The presence of this hard momentum scale allows the 
application of perturbative QCD to predict cross sections for the production of two or more 
high-i?T jets. 



In lowest order (LO) QCD, the hard scattering processes can be classified in two types [jl], 0. 
In the so-called direct process, the photon interacts in a point-like fashion with the quark which 
appears either in the QCD Compton scattering 7g gq or the photon-gluon fusion process 
19 ~^ QQ- the resolved process, the photon acts as a source of partons, which can scatter with 
the partons coming out of the proton. These two processes give two distinctly different types 
of high-ii^T event topologies. The point-like interaction results in events with a characteristic 
(2+l)-jet structure, i.e. two high-E'j' jets and one jet of low Ex fragments of the proton. The 
resolved photon events have a characteristic (2+2)-jet structure, since in addition to the two 
high-ii^T jets and the low-ii^T proton remnant, fragments coming from the photon will produce 
a second \ow-Et jet. While in the (2-|-l)-jet events the total photon energy contributes to the 
hard scattering process, only a fraction participates in the (2+2)-jet events. 

Several methods for detecting the photon remnant jet have been proposed and tested in 
Monte Carlo simulations, generally using a cut on the energy in the photon direction or other 
kinematical variables 0, ^ . Using an Ex cluster algorithm, the ZEUS collaboration was actu- 
ally able to select events with two high-i?T jets and a third cluster in the approximate direction 
of the electron (photon) beam 0. They studied the properties of the photon remnant jet and 
found them to agree with expectations. This shows that photoproduction events with a photon 
remnant jet can be identified experimentally. 



Unfortunately, this does not lead to a separation of the direct and resolved processes, be- 
cause their characterization described above is only valid in LO QCD. In NLO, the direct cross 
section may also have contributions with a photon remnant jet. Therefore its detection does not 
uniquely separate the resolved cross section. In a complete NLO calculation the direct contri- 
butions with a photon remnant jet could be isolated, so that the remaining contributions could 
be compared with the (2-|-l)-jet configuration. Of course, in this case the selection procedure 
in the experimental analysis and in the theoretical calculation must match. To make further 
progress in this direction, a complete NLO calculation of the direct cross section is needed 
which is flexible enough to incorporate the possible selection criteria for the high-i^^ jets and 
the \ow-Et remnant jets. NLO calculations of direct photoproduction have been done in the 
past for inclusive single hadron 0, inclusive single jet 0, H, |^, and for inclusive dijet J^, |l^ 
cross sections. Most of these results are applicable only for special observables in inclusive 
single jet and dijet production. 



The leading order calculation of the direct component involves only two subprocesses: 
10. ~^ 91 and •yg — > qq. By including terms up to O(aaf), new contributions arise. These 
are the one-loop contributions to 'jq gq and •yg qq and the three-body subprocesses 
11 ~^ 199j 11 ~^ nil 11 ~^ ll'l' 1 and yg gqq. Besides the ultraviolet singularities in 
the one-loop contributions, which can be removed by renormalization, additional infrared and 
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collinear singularities occur. They are extracted with the dimensional regularization method. 
It is well known that for suitably defined inclusive observables the infrared singularities in the 
one-loop contributions will cancel against those in the three-body tree diagrams. The collinear 
singularities from initial state radiation and qq production at the photon leg can be factorized 
and absorbed in the respective parton distribution functions. The hard collinear singulari- 
ties from final state radiation will cancel against the corresponding singularities in the one-loop 
graphs. To achieve this technically and to maintain also the necessary flexibility for new choices 
of observables, jet definitions, and experimental cuts, we use the so-called phase space slicing 
method. We introduce an invariant mass cut-off which allows us to separate the phase space 
regions containing the singularities. If y is chosen sufficiently small, the integration over the 
soft and collinear regions can be performed analytically in a straightforward way. Then y is 
just a technical cut which needs not be connected with the experimental jet definition. Of 
course, the suitably defined inclusive one- and two-jet cross sections must be independent of y. 
Our method is very similar to that of Baer et al. [0 who used two distinct cut-off parameters 
for soft and collinear singularities. The phase space slicing method with invariant mass cut-off 
has been used for many processes. The calculation for photoproduction proceeds along the 
same line as the calculation of higher order QCD corrections for multijet cross sections in deep 
inelastic ep scattering first performed by Graudenz ITTf . 



The plan of the paper is as follows: In section 2 we review the LO calculations and fix 
our notation. The techniques to perform the NLO calculations are explained in section 3. 
In section 4 we present some results for inclusive one- and two-jet cross sections in order to 
demonstrate the fiexibility of our method. Section 5 contains a summary of our conclusions. In 
several appendices we provide our analytically calculated results needed for the Monte Carlo 
computations. 



2 Leading Order Cross Sections 
2.1 Photon Spectrum 

The cross sections we have computed are essentially for kinematical conditions as in the HERA 
experiments. There, electrons of Ee = 26.7 GeV produce photons with small virtuality = 
—q^ ~ 0, which then collide with a proton beam of energy Ep = 820 GeV. q = Pe — p'e is 
the momentum transfer of the electron to the photon. The spectrum of the virtual photons is 
approximated by the Weizsacker- Williams formula 

XaF.lAXa) = f (1 + (1 - X^f) In f ^H^~/"^ ) , (1) 

where me is the electron mass and Xa = — E^/ Ee is the fraction of the initial electron energy 
transferred to the photon, p is the four momentum of the incoming proton. In the equivalent 
photon approximation, the cross section for the process ep e'X with arbitrary final state X 
is then given by the convolution 

da{ep e'X) = f dxaF^/e{xa)da{'jp X), (2) 
where dcr(7p — > X) denotes the cross section for jp ^ X with real photons of energy E^ = XaE^. 
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2.2 Jet Cross Sections in LO 



In leading order, the differential cross section for two-jet production takes a very simple form. 
The partonic subprocesses are 7g —>■ gq and 7(7 qq, which we denote generically by 76 — *■ 
PiP2- b is the parton emitted from the proton with momentum po = XbP- Its distribution 
function in the proton Ff,/p{xfy, M^) depends on the momentum fraction Xb and the factorization 
scale Mb. The final state partons have momenta pi and p2, which can be expressed by their 
transverse momentum Et and their rapidities rji and 772- The convention is that the z direction 
is parallel to the proton beam as used in the experimental analysis. From energy and momentum 
conservation one obtains 



Xfi 



^'^ (e^^^+e-"^), (3) 



2E 



Xb = ^(e'^^+e"^)- (4) 

Thus, the kinematical variables of the two jets are related to the scaling variables Xa and Xb- 
For example, the rapidity 772 of the second jet is kinematically fixed by Et, fji, and the photon 
fraction Xa through the relation 

r;2 = -ln(^-e--). (5) 

In the HERA experiments, Xa is restricted to a fixed interval Xa,min < a;^ < max < 1- We shall 
disregard this constraint and allow Xa to vary in the kinematically allowed range Xa^min < a;^ < 1, 
where 

XaMn 2E,Ep-E,ETe^^- ^> 
From equations (3) and (4), we can express Xb as a function of Et, rji, and Xa- 

^ XgEeEre'^' 

2xaEeEp - EpEre-^^ ' ^ ^ 

This function is plotted in Fig. 1 for E^ = 20 GeV and four different values of t]i in the range 
< ?7i < 3 as a function of Xa- Also shown is the envelope of these curves 

'""'''''^ " X EE ' ^^"^ 

which can be obtained from finding the minimum of Xb with respect to rji for fixed Xa- The 
minimal Xa-values for rji = 0, 1, 2, 3 are 0.379, 0.143, 0.0557, and 0.0247. The rapidity of the 
second jet is given by equation (5). 

Depending on Xa, different regions of Xb are explored. Small values of Xb are reached for 
negative rji and the maximal possible Xa, which is Xa = I in our case and Xa = Xa,ma.x in the 
HERA experiments. This is of interest for determining the gluon structure function of the 



proton in the small x region [|12|. Since under experimental conditions Xa is not constrained 
very much, it is possible to achieve minimal values of Xb by selecting the appropriate range in 
Et, rji, and r]2 according to (3) and (4). 
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Figure 1: Contours of constant rji in the Xa — xj, plane for Et — 20 GeV. Xb^env (full curve) is 
the envelope of the curves. 



Of course, studying the gluon structure of the proton this way depends very much on the 
cross sections in the appropriate regions. For direct photoproduction, the two-jet cross section 
for ep — > e' + jet^ + jet2 + X is obtained from 

d^(7 da 

da/dt stands for the differential cross section of the parton subprocess 76 — > piP2- The invari- 
ants of this process are s — (po + qY, t — (po — Pi)^ , and = (po — ^2)^- They can be expressed 
by the final state variables Et, rji, and 772 and the initial state momentum fractions Xa and xi,: 

S = AXaXbEeEp, (10) 

t = -2xaEeETe'^^ = -2xbEpETe-'^\ (11) 
u = -2xaEeETe'^' = -2xiEpETe-'^ . (12) 

So, the dependence of the two-jet cross section on Et, r/i, and 772 is determined through the 
photon distribution function F^/e, the parton distribution functions of the proton F^/p, and the 
cross sections of the hard subprocesses, which depend on s, t, and u. 

For the inclusive one-jet cross section, we must integrate over one of the rapidities in (9). 
We integrate over 772 and transform to the variable Xa using (3). The result is the cross section 
for ep ^ e' -|- jet -|- X, which depends on Et and 77: 



d^a /"^ , ^ , . ^ , ,^2n AE^Et da, , . . ^. 

^ = ? 1,.. dx«-«^Ve(-a)x,F,/,(.„ M,)^^^^^_^^^_ -{lh - p,P2). (13) 



Here, x^ is given by (7) with 771 = 77. 
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The cross sections for the two subprocesses 75 gq and 'jg qq are well known [13 and 
are given by 

^iiqi^gqi) = aasQ^i2CF^ (^-^ - (14) 

-^(7^ QiQi) = aasQi ^ I 7 + ~ ) ' V^^) 

where Cp = {Nq — l)/{2Nc) and Nc is the number of colors. These cross sections are 0{aas). 
The index i denotes the quark flavor and Qi the quark charge. The magnitude of the photon- 
gluon cross section determines whether one is sensitive to the gluon structure function of the 
proton. 



3 Next- To- Leading Order Cross Sections 

The higher order corrections are calculated as usual with the help of dimensional regularization. 
This method regularizes the ultraviolet singularities in the one-loop contributions, which are 
then subtracted using the modified minimal subtraction (MS) scheme. Dimensional regulariza- 
tion can also be used for the infrared and collinear singularities. In the following subsections, 
we shall consider the virtual corrections and the various real corrections needed to obtain finite 
cross sections in the limit dimension n ^ 4. 



3.1 Virtual Corrections up to 0{aal) 

The one-loop diagrams for 75 —>■ gq and 7(7 — > qq have an additional virtual gluon, which leads 
to an extra factor a^. Interfering with the LO diagrams of order 0{aas) produces the virtual 
corrections to the 2 — 2 cross section up to O(aaf). These contributions are well known for 
many years from e^e~ qqg higher order QCD calculations [^4|, |15|. For the corresponding 
photoproduction cross section, one substitutes = and performs the necessary crossings. 



The result can be found in 0, |T^. We have also compared with the results in for deep 
inelastic scattering eq e'gq and eg e'qq, which can be expressed by the invariants s, t, 
and u after setting = 0. The structure of the corrections depends on the color factors Cp, 
CpNci and CpNf, where Nf is the number of flavors in the qq loops. For completeness and for 
later use, we write the final result in the form 

HA... ..0 - - e)g^g (16) 



C^Ms, t, u) - ^NcCfV,2{s, t, u) 



+Cf (^^ + In j QiV; - ^Nc^ T,{s, t, «) | + 0(e), 
Hv{^g^q.q.d = eV/i^^8(l-e)g,^^f^yi^^l^i (17) 



2tt \ s J r(l-2e)4 
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^Nf-^Nc)Tg{s,t,u)\ + 0{e). 
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Hy gives the virtual corrections to the corresponding reactions up to the n-dimensional phase 
space factor 

dt ~ r(l -e)[tu ) 8tts ^ ■' 

and the flux factor l/(2s). The expressions for Vqi, Vq2, Vgi, and Vg2 are collected in appendix 
A. They contain singular terms oc and 1/e (2e = 4 — n), which are always proportional to 
the LO cross sections 

T, ^ (l-.)(-^-i)+2, (19) 

(see eq. (14) and (15)). g is the quark-gluon coupling and fi is the renormalization scale. The 
terms proportional to (|iV/ — y^c) result from the renormalization counter terms. The initial 
gluon spins are averaged by applying a factor of 1/(2(1 — e)). 



3.2 Real Corrections up to 0{aal) 

For the calculation of the hard scattering cross section up to O^aa"^), we must include all dia- 
grams with an additional parton in the final state, i.e. ■jq qgg, •yq — > qqq, 'jq —>■ qq'q', and 
19 The four- vectors of these subprocesses will be labeled by qpo —>■ P1P2P3, where q is 

the momentum of the incoming photon and po is the momentum of the incoming parton. The 
invariants will be denoted by Sij = [pi +Pj)'^ and the previously defined Mandelstam variables 
s, t, and u. For massless partons, the 2 — * 3 contributions contain singularities at Sij = 0. 
They can be extracted with the dimensional regularization method and cancelled against those 
associated with the one-loop contributions. 



To achieve this we go through the following steps. First, the matrix elements for the 2 — 3 
subprocesses are calculated in n dimensions. The corresponding diagrams are shown in Fig. 2 
and are classified as I, II, VII. They are squared and averaged/summed over initial/final 
state spins and colors and can be used for ingoing quarks, antiquarks, and gluons with the 
help of crossing. The diagrams III' account for ghosts to cancel unphysical polarizations of the 
gluons. This classification is the same as in ref. Since the hadronic tensor is symmetric, we 
need to consider only the products IT, IIT, HIT etc. Furthermore, the products VTV, VITV, 
VIITV, VI- V, VII-V, and VII-VI are not singular and vanish with the cut-off describing the 
boundary of the singular region. Such terms will be neglected consistently in the following. 



Next, one can distinguish three classes of singularities in photoproduction {Q^ ~ 0) de- 
pending on which propagators in the squared matrix elements vanish. Examples for these three 
classes are shown in Fig. 3. The X marks the propagator leading to the divergence. In the first 
graph, it is the invariant S23 given by momenta of the final state. Therefore, this singularity 
will be called final state singularity. The second graph becomes singular for Sgi = (q + Pi)'^ = 0, 
when the photon and the final quark momentum are parallel. This is the class of photon initial 
state singularities. In the third graph, the singularity occurs at sqs = (po + PsY, where po is 
the initial parton momentum. This stands for parton initial state singularity. The first class 
is familiar from similar calculations for jet production in e~^e~ annihilation |jT6|, the third class 
from jet production in deep inelastic ep scattering (Q^ 7^ 0) |1T|. The second class occurs only 
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Figure 2: Three-body diagrams for classes I, VII as described in the text. 
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Figure 3: Three-body diagrams with final and initial state singularities. 



for photoproduction [|, |lO|l . 

When squaring the sum of all diagrams in Fig. 2, we encounter terms, where more than one of 
the invariants become singular, e.g. when one of the gluon momenta p3 — so that S23 = and 
So3 = 0. These infrared singularities are disentangled by a partial fractioning decomposition, 
so that every term has only one vanishing denominator. This also allows the separation of the 
different classes of singularities in Fig. 3. Except for additional terms separating the photon 
initial state singularity, we follow the procedure in . It turns out that the results are always 
proportional to the LO cross sections involved in the hard scatter like Tq and Tg in eq. (19) 
and (20) 

Hpjj = KT,. (21) 
Here, F, I, and J denote final state, photon initial state and proton initial state contributions. 



As the last step, the decomposed matrix elements have to be integrated up to Sij < ys. y 
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characterizes the region, where the two partons i and j cannot be resolved. Then, the singular 
kernel K produces terms oc 1/e^ and 1/e, which will cancel against those in the virtual diagrams 
or be absorbed into structure functions, and finite corrections proportional to In^y, Iny, and 
y^. Terms of 0{y) will be neglected. In the following, we shall give the results for the different 
classes of singularities separately. 

3.2.1 Final State Singularities 

In this section, we assume that after partial fractioning the 2 — > 3 matrix elements are singular 
only for Si2 — 0. For the integration over the singular region of phase space, we choose as 
the coordinate system the cm. system of partons pi and p2- The angles of the other parton 
three-momenta po and ps with respect to pi and p2 are shown in Fig. 4. x is the angle between 
Po and p3, 6* is the angle between po and pi, and (p is the azimuthal angle between the planes 
defined by po and pi and po and ps, respectively. 



Po 




Figure 4: Diagram for the three-body final state defining the kinematical variables in the cm. 
system of partons pi and p2- 



In general, the invariants 

s = (Po + g)', (22) 

t = {po-Pi- P2f - 2pi]?2, (23) 

u = (Po-Ps)' (24) 

differ from the corresponding two-body invariants (see eq. (10)-(12)), but approach them for 
P2 = or p2 parallel to pi. Then, t = {po — Pi)^ and u = {po — psY, where pi = pi + P2 is the 
four-momentum of the recombined jet and p^ is the four-momentum of the second jet. This 
produces a jet mass = (pi -|- p2)^ ^ 0, which vanishes in the limit y ^ 0. For p2 — 0, we 
have pi — pi, and for p2 parallel to pi we have pi oc pi. The definitions of t and u are unique 
only in the limit Si2 — > 0. Other choices for the three-parton final state variables differ by finite 
contributions of 0(y), which are neglected throughout. We substitute the angle 9 by 

b=-{l-cose) (25) 
2 
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and introduce the new variable 

= ^. (26) 

Poq 

The three-body phase space in n dimensions can be factorized into two two-body phase 
space terms [|rT 

dPS(') = dPS^'MPS^''), (27) 



where 



and 



with 



dPS(2) 1 /Ans 



dt T(l-e) \tu J Stts 



(28) 



\ s J r(l - 2e) 167r2 I - 2e ^ ^ ^ 



iVf, and A*",^ are normahzation factors: 

N, = /%in-''W- "p'|'/_"^)'' . (32) 

dPS*-^-* is the same phase space factor as in the two-body LO final state. The full range of 
integration in dPS'-''^ is given by z' G [0, — t/s], h e [0, 1], and G [0, vr]. The singular region is 
defined by the requirement that partons pi and p2 are recombined, i.e. p2 = or p2 parallel to 
Pi, so that Si2 = 0. We integrate over this region up to Si2 < ys, which restricts the range of 
z' to < z' < min{— t/ s, y} = yp- 

The matrix elements, which are obtained by squaring the sum of diagrams in Fig. 2, are 
expressed by the variables s, t, u, z', b, and 0. Depending on the structure of the interfering 
diagrams and the color factors, we distinguish seven types of contributions: Fi, F2, F7. Five 
of them are for an incoming quark or antiquark and two for an incoming gluon. The definition 
of the Fi with the corresponding color factors can be read off Tab. 1. 

In class Fi, we have QED-type diagrams proportional to Cp, where a gluon can become soft 
or collinear with a quark. F2 contains QED-type diagrams and interference terms with dia- 
grams containing a three-gluon vertex proportional to CpNc- Class F3 stands for contributions 
of this latter type, where the infrared divergence is caused by the two-gluon final state. F4 has 
only contributions from non-abelian diagrams, and F5 contains the four-quark diagrams. In 
this class, the divergence originates from the collinear qq pair. In Fq and F7, we have contribu- 
tions with incoming gluons proportional to Cp and Nc, respectively. We emphasize that this 
classification rests upon the use of the Feynman gauge for the gluon spin summation with all 
ghost contributions included. Of course, the sum of all terms is gauge invariant. 

After integrating the 2 — > 3 matrix elements with final state singularities ifj?., we obtain 
/dPS<"i/« ^ WS(1 - + 0(.). (33) 
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Class 


Inc. Parton 


Diagrams 


Color Factor 


Fi 
F2 
F3 

F5 
Fe 
F7 


Quark 
Quark 
Quark 
Quark 
Quark 
Gluon 
Gluon 


LI, II-II, II-I 

II I, m i, m il 
mi, mil 
iii-iii 

IV-IV, V-V, VI-VI, VILVII 

I-I, II-II, II-I 
II-I, III-I, III-II 


NcC^p/Nc 
-1/2 NlCp/Nc 
-1/2 NlCp/Nc 

N^Cf/Nc 
1/2 NcCf/Nc 
NcCI/2NcCf 
-1/2 N'^Cf/2NcCf 



Table 1: Classification of the 2^3 matrix elements for final state singularities 

The final results for the Fj are written down in appendix B. They contain the infrared and 
coUinear singularities which cancel against those in the virtual corrections. At this point, it is 
essential that the singular terms are all proportional to the LO matrix elements T^{s,t,u) and 
Tg{s,t,u) and that the variables s, t, and u defined in (22)- (24) correspond to the two-body 
variables in the infrared and coUinear limit with an appropriate definition of the recombined 
jet momentum. 

3.2.2 Photon Initial State Singularities 

To integrate over the singularity Sqs — 0, we use the same pi, p2 cm. coordinate system as in 
the last section. Let us assume that we have isolated the terms in the three-body cross sections 
with a singularity at 5^3 = 0, i.e. where the outgoing quark momentum is coUinear with the 
photon momentum q. In this case, the quark p^ is part of the photon remnant and pi and p2 
are the parton momenta of the two-body subprocess. We introduce a new variable 

z„ = ^e[X„,l], (34) 

where Xa = |^ — Eg/ Eg is now the fraction of the initial electron energy transferred to 
the quark in the hard scatter. Then, Za parametrizes the phase space of the coUinear quark 
P3 = {1 — Za)q and renders the Mandelstam variables in the form 

s = {po + Zaqf, (35) 
t = {Po-Pif, (36) 
u = (P0-P2)'- (37) 



With the additional variable 



the three-body phase space factorizes into 



= ^, (38) 



dPS^^) ^dPS^^MPS^''). (39) 

dPS^^) is identical to (28), and 
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where 



yd4> . T(l-2e) 



and 

H^{z")^{l+^-\ (1-:-—) =l + 0(/') (42) 



can be approximated by 1 as it leads only to negligible terms of 0(|/). The full region of inte- 
gration is given by z" e [0, —u/s], Za G [X^, 1], and e [0, tt], whereas the singular region is 
defined by the requirement that parton is coUinear to g, so that 5^3 < ys/ Za- This restricts 
the range of z" to < z" < min{— ti/s, y} = yj. 

The matrix elements from Fig. 2 can be put into seven classes /i, I2, I7 displayed in Tab. 
2. The classification is again according to incoming parton, color factor, and general structure. 



Class 


Inc. Parton 


Diagrams 


Color Factor 


h 
h 
h 
h 
h 
h 
h 


Quark 

Quark 
Quark 
Quark 
Quark 
Gluon 
Gluon 


I I, n il. II I 
II-I, III-I, III-II 
IV-IV, V-V, VLVI, VILVII 
IV-IV, V-V, VI-VI, VII-VII 
IV IV, V V, VI VI, VII-VII 

I I, n il. II I 

II I, m i, m il 


NcCl/Nc 
-1/2 N^Cf/Nc 
1/2 NcCf/Nc 
1/2 NcCf/Nc 

Cf/Nc 
NcCI/2NcCf 
-1/2 NICf/2NcCf 



Table 2: Classification of the 2 — > 3 matrix elements for photon initial state singularities 

In all classes, the incoming photon can become collinear to a quark line only, since there is 
no direct coupling to a gluon line. Furthermore, there are only single poles, i.e. only collinear 
singularities. The three-body matrix elements are factorized into LO two-body parton-parton 
scattering matrix elements and a singular kernel K. Since the quark (antiquark) produced by 
the photon becomes part of the photon remnant, the remaining two-body process is induced 
by an antiquark (quark) in the initial state, which interacts with a quark or antiquark of the 
same or different flavor or a gluon from the proton. The factorized LO parton-parton scattering 
cross sections are given by the functions [/^(s, u) {i — 1, 2, 3) related to the various processes 
as in Tab. 3. 



The functions Ui{s,t,u) are decomposed with respect to their color structure as follows: 

(43) 



U, - lCF-Nc^){l-e) 



S2 



;i-e)(l.|)-2. 



Nr 

= CFT^a{s,t,u)-^Tu{s,t,u), 
U2 = (44) 




Us = -^il-e){- + e\=^Tsis,t,u). (45) 
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Process 


Matrix element IM]"^ = 8NcCF9*Ui 


m-^ 99 


Ui{s,t,u) 


qq' qq' 


U2is,t,u) 


qq' — > qq' 


U2{u,t, s) 


qq q'q' 


U2{t,s,u) 


qq qq 


U2{s, t, u) + f/2(s, M, t) + Uz{s, t, u) 


qq^ qq 


U2{u, t, s) + U2{u, s, t) + Uz{u, t, s) 


19 19 


-Ui(t,s,u) 



Table 3: LO parton-parton scattering matrix elements 



The functions Tia,Tib,T2, and T3 appear in the final result 



(46) 



which is obtained from the integration of the decomposed matrix elements Hj^ over dPS("\ It IS 
clear that the photon initial state singularities lead to contributions familiar from LO resolved 
photopro duction . 



The terms Jj are collected in appendix C and show explicitly the poles in 1/e proportional 



to 



Pg,^,{Za) = 2NcQ\ 



,1 r 



zt + (1 



(47) 



This function appears in the evolution equation of the photon structure function as an inho- 
mogeneous or so-called point-like term. Therefore, the photon initial state singularities can be 
absorbed into the photon structure function. The necessary steps are well known [jl6|, ITSf. We 
define the renormalized distribution function of a parton a in the electron Fa/e{Xa-, M^) as 



Fa,e{Xa, MI) = / 5,^5(1 - Z,) + — 

J Xa Za L -^vr 
where R has the general form 



Ra^^{Za,M^) 



-^P (z)^^^(^-^X + C iz) 



(48) 
(49) 

jets is then 



with C = in the MS scheme. The renormalized partonic cross section for 76 
calculated from 

dcr(76 jets) = dd-(76 — > jets) - 7^ / dzaRa<-^{za, Ml)da{ah — > jets). (50) 

da{^h jets) is the higher order cross section before the subtraction procedure, and da{ab — >• 
jets) contains the LO parton-parton scattering matrix elements Ui{s,t,u). The factor in 

(49) is combined with the factor in (46) and leads to dependent terms of the form 

'47r/i^V /47ryU^\ 



1 

- -P. 

e 



a<— 7 (•2'a) 



a*— 7 



so that the cross section after subtraction in (50) will depend on the factorization scale 
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(51) 



3.2.3 Parton Initial State Singularities 



The computation of the parton initial state contributions is very similar to the last section. We 
again use the pi, p2 cm. coordinate system as defined in Fig. 4. Singularities occur in this case 
for so3 — 0, when the incoming parton po radiates a soft gluon (pa = 0) or a coUinear parton 
P3. The new variable 

z, = ^e [X,, 1] (52) 

describes the momentum of parton p3 = (1 — Zb)po and the reduced momentum of the parton 
involved in the hard scatter ZbPo- — — Eq/g/Ep is the energy fraction of the incoming 
proton involved in the 2^2 subprocess. The corresponding invariants are 

s = {zbPo + qy, (53) 
t = {zbPo-PiY, (54) 
u = {zbPo-P2f- (55) 



Using the variable 



PoQ 

we can again factorize the three-body phase space into 



= (56) 



dPS(^) = dPS^'MPS^'^), (57) 

where dPS^^^ stands for (28). 

dPSW = f-V ^^^~'\ \ Hb{z"')dfXj, (58) 

\ s J r(l - 2e 167r2 ' ^ ' 



where 



and 



d^, = d."'."'-.^ (^)'^.ta-.^£(l_i4 (59) 
Zb \l-ZbJ r^(l-e) 

Hb{z"') = (1 - z'")-^^^' (l - Y3^) = 1 + 0{z"') (60) 



can again be approximated by 1. The integration over z'" G [0, —u/s], Zb G [Xf,, 1], and (j) G [0, tt] 
is restricted to the singular region of z'" in the range < z'" < min{— m/s, y} = yj. 



The three-body matrix elements can be computed from Fig. 2 and classified according to 

color factor, incoming parton, and type of parton participating in the two-body process as in 
Tab. 4. In the case of final state singularities considered in section 3.2.1, the ingoing parton 
(quark or gluon) was always the ingoing parton of the two-body subprocess. As could already 
been seen in the last section, this is not the case for initial state singularities. For example, if 
an ingoing quark is coUinear to the outgoing gluon (quark), the factorized 2-jet matrix element 
has an initial quark (gluon). 

In Tab. 4, Ji is the class of QED type terms, where a gluon in the final state is collincar 
with the initial state quark. The same occurs in J2, where QED-type and non-abclian diagrams 
interfere. J3 contains the four-quark diagrams, in which an outgoing quark is coUinear with 
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Class 


Inc. Parton 
(3-jet) 


Inc. Parton 
(2-jet) 


Diagrams 


Color Factor 


Jl 
J2 
■h 


Quark 
Quark 

Quark 


Quark 
Quark 
Gluon 


I I, II II, ILI 
II I, III I, m il 

IV-IV, V-V, VI-VI, VILVII 


NcCyNc 
-1/2 NlCp/Nc 
1/2 NcCf/Nc 


Ja 
J5 


Gluon 
Gluon 


Quark 
Quark 


I-I, II-II, II-I 
II-I, III-I, IIMI 


NcCI/2NcCf 
-1/2 NICf/2NcCf 


Jg 


Gluon 


Gluon 


ILL IILL IILII 


-1/2 N^Cf/2NcCf 


Jr 


Gluou 


Gluon 


IILIII 


NICf/2NcCf 



Table 4: Classification of the 2 — > 3 matrix elements for parton initial state singularities 

an ingoing quark. In J4 and J5, we have an ingoing gluon and an outgoing coUinear quark. 
The classes differ only by their color factors. In Jg and J7, an outgoing gluon is collinear with 
the ingoing gluon, where Jq contains QED-type and interference terms and J7 contains pure 
non-abelian terms. 



The next step is to perform the z'" integration with the upper limit yj. The result can be 
written in the form 

,'4.^^Vr(l-«)U_^0(^)_ (61) 



F(l-2e)4 



where the Jj are given in appendix D. Besides s, t, u, and yj, they still depend on the integration 
variable z^. The singular parts take the form 



Jl 
J2 
J3 
Ja 
Je 



[-^NcCF^T,{s,t,u) + 0{e'), 



11 , , , / 1 1 / , -t 3 

--7^^5-.M ^T,{s,t,u) + 0{e'), 

€ Uf J Z 

-'^Pg^giZb) CpTgis, t, u) + 0(e°), 

+ - (--2 + -(^--2-4- C4Nc - U 



eNc 

Nc 



Nc V6 



C|r,(s,i,ii) + 0(e°), (62) 

(63) 
(64) 
(65) 
(66) 



r,(s,t,«) + 0(eO). 



J5 is 0(e) and does not contribute. Jg contains contributions from Jg as well as J7. Some 
of the Jj contain infrared singularities oc 1/e^, which must cancel against the corresponding 
singularities in the virtual contributions. The singular parts are decomposed in such a way 
that the Altarelli-Parisi kernels in four dimensions proportional to 1/e are split off. They also 
appear in the evolution equations for the parton distribution functions and are defined explicitly 
in appendix D. The singular terms proportional to these kernels are absorbed as usual into the 
scale dependent structure functions 



F,/p{X,,M^) 



' dzh 

Xb Zb 



etc 



5,^5(1 -z,) + ^Rl^,,(z,,M^) 



Fb'/p 



(67) 
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where -Fb'/p^^) is the LO structure function before absorption of the colhnear singularities 
and 



2' 



RLA^b, Ml) = --p,M^,) ^^! '\ f ^ 1 + ci^A^,) (68) 



e ' "'r(l -2e) V 



with C" = in the MS scheme. Then, the renormahzed higher order hard scattering cross 
section da{jb — >-jets) is calculated from 

da(76 ^ jets) = da{^b jets) - ^ / dzbR'f,^t,,{zb, M^dai^fb' jets). (69) 

da{'jb — >• jets) is the higher order cross section before the subtraction procedure, and dcr(76' — >• 
jets) contains the lowest order matrix elements Tq{s,t,u) and Tg{s,t,u) in n dimensions. This 



well known factorization prescription |]T8| removes finally all remaining collinear singularities. 
It is universal and leads for all processes to the same definition of structure functions if the 
choice concerning the regular function C in (68) is kept fixed. Similar to the case of photon 
initial state singularities, the higher order cross sections in (69) will depend on the factorization 

scale Mfe due to terms of the form Ph^i,/{zi,) In 



3.3 Jet Cross Sections in NLO 

To obtain a finite cross section, we must add the four parts considered in sections 3.1 and 
3.2.1-3.2.3. Then, the poles in 1/e must cancel, and we can take the limit e — 0. The result is a 
special kind of two-jet cross section, where the recombination of two partons into one jet or the 
recombination of a parton with the photon or proton remnant jet is done with an invariant mass 
cut-off y. In LO, this cross section was written down in (9). Including the NLO corrections, it 
has the form 



dV 



dE^dr]idr]2 



o.da 

+XbFb/p{xb, M^) — [ab P1P2) 



da , 



da 
dt 



(7^ 



P1P2) 
■ P1P2) 



(70) 



In (70), ^{iQ ~^ P1P2) and ^{l9 ^ P1P2) stand for two-body contributions in LO and NLO 
together with analytically integrated contributions of the soft and collinear divergent regions 
of the three-parton final state. The contributions from the photon initial state singularities are 

denoted ^{o-b — >• piP2)- ^il^ ~^ P1P2) and ^{l9 P1P2) originate from the integration of 
the soft and collinear singularities of the parton initial state. 



The quark induced cross section can be written as 



^{iq - P1P2) = a,{l^')CT, + ^Mlc{T,Ag + B, 



(71) 



where 



1 



s N 



A = --(lliVc-2iV^)ln-^ + ^(lny^-^) -iVc^lny^ 



q g V — o --/^- 2 ■ 3 gy - - e 
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(72) 
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+ {2Cf 
+ Nc 

+ Cf 



Nc) 



- In^ - - ^ In^ — + In yp In — 



4 

yps 



-t 

s 



Lie 



1, 2 1 

m Vp ___ 

2 2 -u 



+ - In^ Li 

4 s 



yFS\ 67 TT^ 
^ ^ ~J + 18 " T 



TT^ 7 2 3 



T = 



Q 11 1 

_3C^_ln (2C^_iVc,) 

It s 2 

+ f2 + ^)ln^^ 
\ UJ s 

, and C = -^QfCf 

s u 



2 In 



—MS 



u 



For the gluon induced cross section, we have 

3 



da 



-C{T,A, + B,), 



where 

An 



In^ — + In^ — - 2 In^ - 3 Inyp 



^{nNc-2Nf) In^ + Cp 

^ 1 2 1 1 2 2 , 1 1 2 

In — ^ m h In yp + - in 



(73) 

(74) 

(75) 
(76) 



1 , 9 

-In — 



6 



+Nc 



-t 



-t 1 

- + T 



In^ 



-U 



B„ — 3Ce 



t . —u u . —t\ 1 

_ln + ln_ + - {2Cf 

u s t s J 2 



Nc) 



' , tU ( ^ U\ , n —U 



T 



+ (2 + -)ln2 — 

t ^u 

u t 



(77) 



(78) 



The cut dependence of the final state NLO corrections in Ag and Ag is contained in the yp 
dependent terms. For yp — 0, these terms behave hke {— In' yp), which leads to unphysical 
negative cross sections for very small yp. Thus if y is used as a physical cut, it must be suf- 
ficiently large. In most of the applications, we shall use these results for computing inclusive 
cross sections, in which the y dependence of the two-jet cross section cancels against the y de- 
pendence of the numerically calculated three-jet cross section. The quadratic terms in (Iny^) 
are connected with the soft divergence of the three-parton final state. 



The contributions from the photon initial state singularities necessarily have a photon rem- 
nant. They depend on the cut-off yj only through terms proportional to (Inyj). This contri- 
bution can be considered separately only if yj is a physical cut, i.e. if it is sufficiently large. 
In case one introduces other kinematical conditions to define the photon remnant terms in the 
NLO direct cross section, one must add the appropriate three-body contributions. The two-jet 
cross section is calculated from 



dt ^ 



1 dZa 



f 

Jx, 



a 

Pqi^liZa) In 



(79) 



Za)yis' 



ZaMl 



l\+2Nc 
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In (79), Ma is the factorization scale at the photon leg. This dependence must cancel against 
the Ma dependence of the LO resolved contribution (see eq. (48)). This has been demon- 
strated explicitly in for the inclusive single jet cross section. The ^{o-b — > P1P2) are the 
cross sections for the 2 —>■ 2 subprocesses which were given in section 3.2.2 in terms of the 
functions Ui{s,t,u). These are the well known LO cross sections also used for the calculation 
of hadron-hadron cross sections. 



From the parton initial state singularities, we deduce 



da 
dt 



iiq P1P2) 



da 
dt 



iig P1P2) 



XbFqIp 



XbFg/p 




— — XbFqJp 



Zb 

Xb_ 

Zb ' 



Ml 



CpJliZb) + Pq^q{Zb) In 



Pq^g{Zb) I In 



Zb)yjs^ 



ZbM, 



Nc . 



MZb)+ Pg^g{Zb)\Ylj^ 



(80) 



-1 + 




(81) 



The functions Jj are the regular parts of the functions Ji defined in appendix D. They contain 
the cut-off yj and are distributions in Zb. The dominant terms for ?/j — are again propor- 
tional to (In^ yj) and are related to the soft divergence in the parton initial state. J2 does not 
contribute after the appropriate crossing of {t u). 



The results in (70)-(81) are equivalent to those of Baer et al. [0. They cannot be compared 
since these authors introduce a different cut for integrating over the soft singularities. Instead 
of invariant masses, they use a cut on the energy of the gluon. Terms that are independent of 
soft singularities agree. For example, this is the case for the cross section in (79) coming from 
the collinear photon initial state singularity. 



4 Inclusive One- and Two-Jet Cross Sections 

In this section, we present some characteristic numerical results for one- and two-jet cross sec- 
tions which have been obtained with our method of slicing the phase space with invariant mass 
cuts. In a short communication we have used this method to calculate the differential two-jet 
cross section as a function of average rapidity fj = ^{rji -1-772) of the two jets of highest transverse 
energy with the cut on \rii — ri2\ < 0.5 and compared it to recent experimental data of the 
ZEUS collaboration at HERA [^. For this paper, we have calculated various one- and two-jet 



distributions without applying special cuts on kinematical variables of the initial or final states 
dictated by the experimental analysis, although our approach is particularly suitable for this 
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as has been shown in 



20 



The input for our calculations is listed in the following. We have chosen the CTEQ3M 



proton structure function |2l|, which is a NLO parametrization with MS factorization and 
A*^"'^ = 239 MeV. This A value is also used to calculate the two-loop as value at the scale 
H = Ej-. The factorization scale is also M = E^- We do not apply a special cut on the photon 
energy fraction Xa in (2) but integrate from Xa,min to 1. Q^ax — ^ GeV^ is the same as used in 



the jet analysis of the ZEUS experiment [22 



The further calculation is based on two separate contributions - a set of two-body contri- 
butions and a set of three-body contributions. Each set is completely finite, as all singularities 
have been cancelled or absorbed into structure functions. Each part depends separately on the 
cut-off y. If y is chosen large enough, the two parts determine physically well defined two-jet 
and three-jet cross sections. Our analytic calculations are valid only for very small y, since 
terms 0{y) have been neglected in the analytic integrations. For very small y, the two pieces 
have no physical meaning. In this case, the (Iny) terms force the two-body contributions to 
become negative, whereas the three-body cross sections are large and positive. In pO|, we have 
plotted such cross sections for y = 10~^. However, when the two- and three-body contributions 
are superimposed to yield a suitable inclusive cross section, as for example the inclusive single- 
jet cross section, the dependence on the cut-off y will cancel. Then, the separation of the two 
contributions with the cut-off y is only a technical device. The cut-off only serves to distinguish 
the phase space regions, where the integrations are done analytically, from those where they 
are done numerically. Furthermore, y must be chosen sufficiently small so that experimental 
cuts imposed on kinematical variables do not interfere with the cancellation of the y dependence. 

First, we consider the inclusive single-jet cross section. To achieve this, we must choose 
a jet definition, which recombines two nearly collinear partons. As usual, we adopt the jet 
definition of the Snowmass meeting . According to this definition, two partons i and j are 

recombined if Rij < R, where Ri = \J (r]i — r]jy + {(pi — (pjY and ?7j, 0j are the rapidity and 
the azimuthal angle of the recombined jet. We choose R = 1 m all of the following results. 
This means that two partons are considered as two separate jets or as a single jet depending 
whether they lie outside or inside the cone with radius R around the jet momentum. In some 
cases, it may occur that two partons i and j qualify both as two individual jets i and j and as 



a combined jet ij. In this case, we make no further selection as it is done in p4|, where only 
the combined jet is counted. In NLO, the final state may consist of two or three jets. The 
three-jet sample consists of all three-body contributions, which do not fulfill the cone condition. 

In Fig. 5, the results for the inclusive one-jet cross section d^cr/d^^rdr^ are shown for various 
rapidities = 0, 1, and 2 as a function of Et- The NLO and the LO predictions are plotted. 
The LO curve is calculated with the same proton structure function and the same as the 
NLO curve. Only the hard scattering cross section is calculated in LO. For i? = 1, these two 
cross sections differ only very little. Of course, the NLO cross section depends on i?, whereas 
the LO curve does not. To separate the curves for the three rj values, the result for 77 = (?7 = 2) 
has been multiplied by 0.1 (0.5). 

The maximum of the one-jet cross section is near r] = 1. For Et = 20 GeV, this can be seen 
from Fig. 6, where d a/dExdr] is plotted as a function of i] between 77 = —1 and r] ~ 4.2. At this 
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ri = 0, LO 
77 = 0, NLO 
77 = 1, LO 
77 = 1, NLO 
77 = 2, LO 
7j = 2, NLO 
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X 0.1 
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Figure 5: Inclusive single-jet cross section d^a/dExdr] as a function of Et for various rapidities 
?7 = 0, 1, 2 in LO and NLO. The cross section for 77 = (77 = 2) is multiplied by 0.1 (0.5). 
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value of Et, the kinematical boundaries, where the cross section goes to zero, differ only little in 
LO and NLO. The /c-factor for the NLO prediction is A; ~ 1.2 for the central region of < 77 < 3. 

Our single-jet distributions have been compared with the earlier results of Bodeker 0, who 
used the subtraction method to cancel soft and coUinear singularities. We found perfect agree- 
ment. To achieve this, we had to use a rather small value of y = 10^^. Furthermore, we had 
to add some finite contributions of 0(1) for y —>■ 0, which correspond to two-jet contributions 
where an outgoing parton is coUinear to another outgoing or an ingoing parton without the 
corresponding pole term in the 2 — > 3 matrix elements after partial fractioning. We also checked 
that our results are independent of y, if y is sufficiently small. 

The more essential results are on inclusive two-jet cross sections shown in the following fig- 
ures. In Fig. 7, we present d^a/dE'j'^driidri2 as a function of Et^ for rji = 1 and various choices 
of r]2 = 0, 1, 2. Here, Et^ and rji are the transverse energy and the rapidity of the so-called 
trigger jet. 772 is the rapidity of the second jet, so that Et^ and Et2 are the jets with the highest 
transverse momenta for the three-jet contribution. For exactly two jets in the final state, we 
have Eti = Et^- In Fig. 7, we can see how the cross section decreases when 772 is chosen away 
from the maximum region at 772 = 1- ?7i = 1 is always kept fixed. To disentangle the three 
curves, we have multiplied the 772 = 2 (772 = 0) cross section by a factor of 0.5 (0.1) as in Fig. 5. 
The LO and NLO predictions are very much parallel. The fc-factor is A; ~ 1.2, similar as in the 
one-jet cross section. The only large difference between the LO and NLO cross sections occurs 
at the boundary of the phase space. For example, the cross section for 772 = has an additional 
tail in the Et^ distribution coming from the three-jet cross section, where all three partons lie 
outside the cone with radius R = 1. Such plots can be produced also for other choices of r/i 
outside the region where the cross section is maximal. 

We have studied the inclusive two-jet cross section also as a function of rji and 772 for fixed 
Et^. As an example, we show in Fig. 8 the two-dimensional distribution d? a / dET^driidri2 for 
Et^ = 20 GeV in form of a lego-plot in the intervals 771, 772 G [—1, 4]. Over the whole range of rji 
and 7/2, the cross section is approximately symmetric in rji and 7/2. We show this cross section 
in LO and NLO, and we see clear differences between the LO and the NLO prediction. 

This becomes even clearer when we plot the projections for fixed t]i as a function of 772. In 
Fig. 9, these projections are shown for 7/1 = 0, 1, 2, and 3 for LO and NLO. We see that in 
NLO the cross section is always larger than the LO cross section which is consistent with the 
A;-factor larger than 1 given above. 

It is clear that many more distributions or partially integrated cross sections using other 
two-jet variables can be calculated with the phase space slicing method. As a final application, 
we have studied the R and the y dependence of the two-jet cross section, integrated over 771,7^2 
and Et^ with a minimal -ETi.min = 10 GeV. The result is shown in Fig. 10, where the LO 
cross section is independent of R. The NLO cross section is the sum of positive three-body 
contributions and negative two-body contributions, both calculated with the invariant mass 
cut y = 10~^, and this is not independent of R. However, the dependence is rather weak. The 
cross section increases with increasing R from i? = 0.1 to i? = 1 by a factor of 2 and lies above 
the LO result for R > 0.6. For R < 0.1, the y cut becomes effective. Therefore we have plotted 
the curves only for R > 0.1. The sum of the two-body and three-body contributions depends 
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Figure 7: Inclusive dijet cross section d^a/dETidr]id7]2 as a function of Et^ for ?7i = 1 and three 
values of 772 = 0, 1, 2. The cross section for 772 = (772 = 2) is multiphed by 0.1 (0.5). 




Figure 8: The LO and NLO triple differential inclusive dijet cross section for E^^ = 20 GeV as 
a function of rji and 772- 
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7;'i = 3, LO 

7?, = 3, NLO 




Figure 9: Projections of the LO and NLO triple differential inclusive dijet cross section for 
= 20 GeV as a function of 772 for fixed rji = 0, 1, 2, and 3. 



on R like a + bin R + cR^ due to the recombination of two parton momenta Pi and pj inside 
the cone with radius R. This way, the singled out transverse energy of a jet can differ from the 
individual parton E^s. This R dependence has been studied before for the inclusive one-jet 
cross section with similar results p. 

We have also calculated the same cross section as a function of the invariant mass cut y 
without the cone recombination. In this case, the cut y operates as a physical cut. Of course, 
the results for the two-body contributions (being negative for y < 3 ■ 10^^) and the three-body 
contributions (being above the sum for y < 3 - 10~^) only make sense for sufficiently large values 
of the y cut. For ?/ — > 0, the two-body (three-body) contributions decrease (increase) like — In^ y 
(In^y). The results are presented in Fig. 11. Here we see that the sum of both contributions 
depends only moderately on y (oc Iny) due to the recombination of partons in the final state. 
This dependence on y of the summed contribution is equivalent to the R dependence of the in- 
tegrated cross section shown in Fig. 10. We remark that in Fig. 11 the respective cross sections 
must still be corrected by finite terms of 0{y). The two-body contribution is calculated from 
the analytical formulas given above, where these terms are neglected. This should influence 
our results only for y > 2- 10~^. It is, however, no problem to calculate these terms numerically. 
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Figure 10: Integrated inclusive dijet cross section for i?Ti,min = 10 GeV as a function of R in 
NLO (dashed-dotted) and LO (full). The 2 body (dashed) and 3 body (dotted) contributions 
with y = 10^^ are plotted separately. 
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Figure 11: Integrated inclusive dijet cross section for -BTi.min = 10 GeV as a function of y in 
NLO (dashed-dotted) and LO (full). The 2 body (dashed) and 3 body (dotted) contributions 
are plotted separately. 
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5 Summary 



Differential cross sections d^a/dE^driidri2 have been calculated in NLO for direct photopro- 
duction. Infrared and coUinear singularities are cancelled with the phase space slicing method 
using an invariant mass cut-off. This method allows to incorporate various cuts on the final 
state dictated by the analysis of experimental data and to perform calculations for different 
choices of jet algorithms. Analytical formulas for the different contributions giving the depen- 
dence on the slicing parameter are derived. This is needed if one wants to develop a complete 
description of photoproduction events including hadronization and other corrections. Numeri- 
cal results for the two-jet inclusive cross sections at HERA have been presented, and LO and 
NLO predictions have been compared. For a cone radius of R = 1, the NLO corrections lead 
to an increase of the order of 20% compared to the LO prediction which is calculated with the 
same proton structure function and the same as the NLO prediction. 
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A Virtual Corrections 



In this appendix, we give the interference terms of the virtual corrections with the LO QCD 
Compton scattering •yq — gq (= Tq{s,t,u)) and the LO photon-gluon fusion ■jg qq (= 
Tg{s,t,u)) matrix elements as defined in equations (16) and (17). The results depend only on 
the Mandelstam variables s, t, and u: 
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B Final State Corrections 

This appendix contains the real corrections coming from the integration of the 2 — > 3 matrix 
elements over the final state singular region of phase space. In addition to the Mandelstam 
variables, the result depends on the integration cut-off yp'- 
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C Photon Initial State Corrections 



Like gluons, real photons can split into qq pairs. After integration over the collinear region of 
phase space, the Born matrix elements for qq — > gg (= Tia^b{s , t , u)) , qq' qq' (= T2{s,t,u)), 
and qq qq (= T2^3{s,t,uj) can be factorized out. The result depends on the Mandelstam 
variables, on the integration cut-off yi, and on the integration variable Za'- 
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The Altarelli-Parisi splitting function Pg^^-y is equal to 2NcQlPq^g (see Appendix D). 
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D Parton Initial State Corrections 



In this appendix, the parton initial state corrections are given as a function of the Mandelstam 
variables s, t, u, of the cut-off yj, and of the additional integration variable Zb. 
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for convenience. The AltareUi-Parisi sphtting functions are defined as 
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As the integration over runs from X;, to 1, the + distributions in this paper are defined as 



DAg] 



where 
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(D.13) 



and h{zb) is a regular function of Zb This leads to additional terms not given here explicitly 
when (D.12) is transformed so that both integrals are calculated in the range [X^, 1]. 
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